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5 ■ ABSTRACT 

'NT I 

04 I We investigate the transfer of power between different scales and coupling of 

■ modes during non-linear evolution of gravitational clustering in an expand- 
> ■ 

' ing universe. We start with a power spectrum of density fluctuations that 



o 



is exponentially damped outside a narrow range of scales and use numerical 



I simulations to study evolution of this power spectrum. Non-Linear effects gen- 

. erate power at other scales with most power flowing from larger to smaller 

' scales. The "cascade" of power leads to equipartition of energy at smaller 

Oh' scales, implying a power spectrum with index n « — 1. We find that such a 

O I spectrum is produced in the range 1 < S < 200 for density contrast 6. This 

C/3 . result continues to hold even when small scale power is added to the initial 

power spectrum. Semi-analytic models for gravitational clustering suggest a 
' tendency for the effective index to move towards a critical index rtc ~ —1 in 

^ I this range. For n < ric, power in this range grows faster than linear rate, while 

if n > Uc, it grows at a slower rate - thereby changing the index closer to ric. 
At scales larger than the narrow range of scales with initial power, a k'^ tail 
is produced. We demonstrate that non-linear small scales do not effect the 
growth of perturbations at larger scales. 
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2 J.S.Bagla and T.Padmanabhan 
1 INTRODUCTION 

Gravitational instability leads to growth of density inhomogeneities in an expanding uni- 
verse. In Fourier space, one can study this growth as the evolution of Fourier modes of 
density contrast, which evolve independently of each other in the linear regime. The actual 
growth rate depends on the background cosmology and in a matter dominated universe with 
r2 = 1, the amplitude of Fourier modes grows as the expansion factor a{t) oc t^/s [Th is is the 
background model considered in this paper]. 

The situation is quite different if the linear perturbation theory is not applicable and 
one has to consider the effect of coupling between different modes. The key effect of such a 
coupling will be the transfer of power between different length scales. Some aspects of power 
transfer have been studied by using N-Body simulations in the past. For example. Little, 
Weinberg and Park ( p.991|) conducted a series of N-Body experiments with an initial power 
spectrum that was truncated at different length scales. They concluded that the structure 
and appearance of a non-linear universe is dominated by a small range of scales, centered 
around the scale that is becoming non-linear. A detailed study of transfer of power and 
evolution of phases was carried out by Soda and Suto (|1992| ) using Zeldovich approximation, 
wherein they showed that one dimensional collapse leads to equipartition of fluctuation power 
per unit wavenumber. They also carried out a few experiments with N-Body simulations for 
gravitational clustering in three dimensions and noted that a similar trend is indicated in 



this well. Klypin and Melott ( |1992| ) studied evolution of ratio of kinetic energy 



at different scales using N-Body simulations and concluded that this ratio tends towards a 
universal value in non-linear evolution of clustering and approaches the value for the n = — 1 
power spectrum. 

The simplest context in which one can view power transfer is the following: Some amount 
of power is injected at a given scale at t = ti and the nonlinear coupling is allowed to 
transfer it to other scales. By studying such an evolution numerically and computing the 
power spectra at later epochs we will be able to understand how the transfer of power takes 
place. In this paper, we will demonstrate a few generic features of power transfer between 
different modes by studying evolution of power spectrum in some such models and present 
theoretical arguments regarding their validity. The theoretical arguments are based on some 
of the recent approaches to the study of nonlinear clustering [briefly discussed in the next 
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section] and are presented in the spirit of a paradigm to understand the numerical results. 
They should, of course, not be thought of as rigorous mathematical proofs. 



2 TRANSFER OF POWER 

The evolution of the density contrast in Fourier space can be described by an equation of 
the form ([Peebles 1974| ) 



+ 2-5k = 47rGp5k + Q (1) 
a 

where 6k{t) is the Fourier transform of the density contrast, p is the background density 
and Q is a non-local, non-linear function which couples the mode k to all other modes k'. 
Coupling between different modes is significant in two cases: (i) An obvious case is one 
with > 1, i.e. the amplitude of density perturbations at the scale of interest is of order 
unity or larger, (ii) A more interesting possibility arises for modes with no initial power [or 
exponentially small power]. In this case nonlinear coupling provides the only driving terms, 
represented by Q in eqn.(|I]). These generate power at the scale k through mode-coupling, 
provided initial power exists at some other scale. Note that the growth of power at the scale 
k will now be governed purely by nonlinear effects even though 5^ ^ 1 ■ As we shall see, this 
fact leads to some interesting effects. 

The exact solution to eqn.(|l]) is, of course, not known. But it is possible to understand 
some simple features of nonlinear clustering by studying the characteristics of partial dif- 
ferential equations that govern the evolution of two-point correlation function. [See Peebles 
( |1980|) and Nityananda and Padmanabhan ( [L994| ).] One may draw the following conclusions 



from such a study for hierarchical models [For more details see Padmanabhan ( p.996a|) and 
( I1996BD .] : 



• The power at a scale x, at an epoch a is related closely to the linearly extrapolated 
power, at a scale / where 

l = x{l + e(a, . ^^^^ x) = ^r e(a, y)y''dy (2) 



Jo 



Here ^ is the two point correlation function and a is the scale factor. ^ (x, a) is the mean 
correlation averaged up to the scale x at the epoch a ( Hamilton et al 1991| ). 



• N-Body simulations of hierarchical models show that the non-linear mean correlation 
function at x relates to the linear mean correlation function at / in an almost universal 
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manner ( |Hamilton et al 1991| ). The relation can be approximated by the following map 
( Pagla and Padmanabhan 1995|) 

a(a,/) (a<1.2, e<1.2) 

0.7^L(a, If (1.2 < < 6.5, 1.2 < f < 195) (3) 



11.7a(a,0'^' (6.5<a, 195<0 



provided we can assume stable clustering of virialised objects in the non-linear regime. 

Relation between ^(a, x) and originally given by Hamilton et al. ( 1991| ) based 
on the N-body simulation data. The eqn.(^) is an approximate fit to N-Body data which 
has the advantage of being a piecewise power law fit in the three regimes. The three regimes 
used in this power law fit have some physical relevance as has been shown in a recent paper 
( padmanabhan 1996a| ) that also justifies the slopes of the power laws used in each regime. 
Several authors have provided more exact fitting functions to describe the three regimes in a 
unified manner. [See Hamilton et al. ( |1991|) ; Peacock and Dodds ( |1996|) ; Jain, Mo and White 
( |1995|) .] There is also some controversy regarding the actual index in the intermediate regime 
^ ^i] ^"^^ whether the relation given above is truly "universal" [see for example Peacock 
and Dodds ( |1996|) ; Jain, Mo and White (p95D ; Padmanabhan et al. (p96|) .]. However, all 
the models suggested in the literature lead to the following feature : In the quasi-linear 
regime, there exists a critical index such that the growth of power is faster than for 
spectra with n < ric and slower than a? for spectra with n > ric- 

With the simple scaling given above, it is easy to show that Uc = —1; but if a more 
accurate fitting function is used then this value may vary around —1. For the sake of illus- 
tration we shall take rif. = —1 but the general arguments in this paper do not depend on the 
particular choice of ric- The existence of such an index suggests that, during the evolution, 
there will be a tendency for the power spectra to acquire such an index. This conclusion - 
which is a direct consequence of the results described above - is worth testing in numerical 
simulations. 



3 NUMERICAL EXPERIMENTS 

We shall now try to see what these results might imply for the transfer of power in a more 
general context. But before describing the results of such an experiment, we would like to 
spell out the theoretical expectations, this is important because it shows - up front - what 
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the results should be and helps us to understand them. Needless to say, the real justification 
for the claims made below comes from the results of the simulations. 

To begin with, it is well known that the power transfer in gravitational clustering is 
mostly from large scales to small scales. [A significant exception is the generation of the /c^ 
tail which we shall discuss towards the end of the paper.] This is clearly borne out by the 
numerical experiments of Little, Wienberg and Park ( |1991|) that showed that the structure 
and appearance of nonlinear structures in simulations is largely independent of the initial 
power at small scales. Suppose we start with a power spectrum that is centered at some 
scale ko = 2tt/Lq and has a small width Ak. First structures to form in such a system are 
voids with a typical diameter Lq. Formation and fragmentation of sheets bounding the voids 
leads to generation of power at scales L < Lq. First bound structures form at the mass scale 
corresponding to Lq. In such a model the linear ,^ at L < Lq is nearly constant with an 
effective index of n —3. Assuming we can use eqn.(^ with the local index in this case, 
we expect the power to grow very rapidly as compared to the linear rate of a^. [The rate 
of growth is a® for n = —3 and for n = —2.5.] Different rate of growth for regions with 
different local index will lead to steepening of the power spectrum with an accompanying 
slowing down of the rate of growth. This rapid growth is expected in the quasi-linear regime 
and should lead to a power spectrum with the critical index as its slope. 

Consider next a more complex situation, with initial power concentrated around two 
scales Lq and Li < Lq. If we assume that the power at Li has a higher amplitude then 
the smaller scale(s) will reach the quasi-linear phase before the larger one(s) and - in the 
subsequent evolution - will approach the critical index n^. We again expect the spectrum 
to have n = Uc a.t scale L < Li. If the largest non-linear scales dominate the non-linear 
evolution of perturbations, as shown in a limited context ( [Little, Weinberg and Park 19911) , 
then the spectrum at scales Li < L < Lq should also approach one with an index n = 
after the larger scale becomes nonlinear. However, if the extra small scale power implies a 
much steeper spectrum, we should expect the small scale power to grow at a slower rate in 
order to match up with the n = Uc power spectrum at larger scales. If this is indeed the 
case, then the power spectrum is effectively driven by the largest scale which is entering 
the quasi-linear phase at the epoch of consideration. This scale is expected to influence the 
scales in the quasi-linear regime. 

The question we would like to address in this paper is : To what extent is the above 
qualitative picture supported by numerical simulations ? To answer this question we will 
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6 J.S.Bagla and T.Padmanabhan 

use N-Body simulations and study evolution of some toy models. We begin with a brief 
summary of the N-Body code used for these numerical experiments. For greater details of 
the PM code used here see Bagla and Padmanabhan (1996) 

All simulations used a PM [Particle Mesh] code and (128)^ particles in a (128)3 box. In 
the units of length used here, each side of the simulation box measures 128 units. We used 
the TSC [Triangular Shaped Cloud] for interpolation and the "poor man's" Poisson solver 
for solving the Poisson equation in Fourier space. Force was computed in Fourier space from 
the potential. For more details on PM codes, see Hockney and Eastwood (|1980| ). Numerical 
integration of the equation of motion was done using the standard Leap-Frog method. Step 
size for integration was chosen by enforcing an upper limit on the maximum displacement of 
a particle in one step. The "woe factor" [lOAt/t^yn] for simulations used here is about 0.2. 

In all the figures for power spectrum in this paper, we shall show only that region in the 
/c-space where the uncertainty is small. The smallest scale shown in these figures is 1.5Lny 
where Lny is the Nyquist scale. Average number of particles enclosed in a cube of this size, or 
larger, contains sufficiently large number of particles and hence the error in power spectrum 
is fairly small at these scales. The largest scale used in these plots is L = L^ox/'^', thus 
the averaging for power spectrum is done over at least eight independent regions in the 
simulation box. 

We used three models for our study. Parameters of these models were chosen as follows 

• Initial power spectrum for model I, the "reference" model, was a Gaussian peaked at 
the scale ko = 2tt/Lq]Lo = 24 and having a spread Ak = 27r/128. The amplitude of the 
peak was chosen so that A;j„(A;o = 27r/Lo,a = 0.25) = 1, where A'^{k) = k^P{k) / (2'k'^) and 
P{k) is the power spectrum. Needless to say, the simulation starts while the peak of the 
Gaussian is in the linear regime (A(A;o) *C 1). 

• Model II had initial power concentrated in two narrow windows in fc-space. In addition 
to power around Lq = 24 as in model I, we added power at ki = 2tt/Li;Li = 8 using a 
Gaussian with same width as that used in model I. Amplitude at Li was chosen five times 
higher than that at Lq = 24, thus A;j„(fci, a = 0.05) = 1. 

• Model III was similar to model II, with the small scale peak shifted to ki = 271 / Li; Li = 
12. The amplitude of the small scale peak was the same as in Model II. 

We now describe results of simulations of these models. Panels of figure 1 show a slice 
from the simulation volume at four epochs for each model. A comparison of non-linear 
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Figure 1. This figure shows a shce from simulations of model I at four epochs. The upper left panel corresponds to a = 0.25, 
the upper right panel to a = 0.5, the lower left panel to a = 1 and the lower right panel to a = 2. The thickness of the slice is 
6L and the width and height are 128L where L is one grid length. 

structures in these panels shows that at early epochs the appearance of non-linear objects 
in these models is very different from each other. Dissimilarities between different models 
slowly disappear and the appearance at later epochs is dominated by the larger wave-mode. 
This is true even for model III that has a lot of small scale power. [If we define an effective 
index of the power spectrum by joining the peaks of two Gaussian-s then the effective index 
Ueff ~ for model II and n^jf ~ 2 for model III.] It is clear from these pictures that the 
large scale mode takes a longer time to assert itself for the case with greater small scale 
power. 

To compare the evolution of these models in a more quantitative fashion we use the 
power spectrum. Top panel of figure 2 shows evolution of power spectrum for model I. The 
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Figure 1. Continued. This figure shows the corresponding sHce from simulations of model II. 

y-axis is A{k)/a, the power per logarithmic scale divided by the linear growth factor. This 
is plotted as a function of scale L = 211/ k for different values of scale factor a{t), curves 
are labeled by the value of a. As we have divided the power spectrum by its linear rate 
of growth, the change of shape of the spectrum occurs strictly because of non-linear mode 
coupling. It is clear from this figure that power at small scales grows rapidly and saturates 
to growth at a rate close to the linear rate [shown by crowding of curves] at later epochs. 
The effective index for the power spectrum approaches n = — 1 within the accuracy of the 
simulations. Thus this figure clearly demonstrates the general features we expected from our 
understanding of scaling relations. 

The other two panels of figure 2 show the corresponding curves for models II and III 
respectively. These models had power at small scales in addition to the power at large scales. 
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Figure 1. Continued. This figure sliows the corresponding sHce from simulations of model III. 

The large scale power is same in all the models and the initial conditions for the relevant 
region in fc-space had same initial phases, making comparison meaningful. These figures 
show that when large scales become non-linear, power at small scales grows at a rate slower 
than the linear rate in the quasi-linear regime, till we get a power spectrum with n ^ —1 
for L < Lq. The amplitude of power spectrum for models II and III at this stage is same as 
the corresponding power spectrum in model I. A comparison of the lower panels also shows 
that the approach to Uc is slower for the model with more small scale power (model III). 

Figure 3 shows power spectra of all three models at a late epoch. At this epoch A;j„(A;o) = 
4.5 and it is clear from this figure that the power spectra of these models are very similar 
to one another. 

The growth of power at three different scales, L = 8, 12, 24, is plotted for the three models 
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L=27T/k 



Figure 2. This figure shows evolution of power spectra for the three models. Thick lines show the linear power spectrum for 
these models and thin lines show the non-linear power spectrum from N-Body simulations. The y-axis is square root of power 
per logarithmic scale divided by a. With linear growth rate divided out, only non-linear evolution can modify the spectrum in 
this plot. X-axis is the length scale. In model I (top panel) the initial power spectrum is a Gaussian peaked at a length scale 
of Lo = 24 and amplitude adjusted to make it reach non-linearity at a = 0.25. Power spectra for different epochs are labeled 
by the scale factor. This figure demonstrates that power is generated at smaller scales and at late times this power saturates 
to a power spectrum with n ^ —1. The second panel shows the evolution of model II. In this model there is additional power 
(apart from what is there in model I) at Li = 8 with an amplitude five times higher. 
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L=27T/k 

Figure 2. Continued. This panel shows the same curves for model III. In this model there is additional power (apart from 
what is there in model I) at Li = 12 with an amplitude five times higher. Notice that at late times power spectra of all three 
models are very similar. 



O 




5 10 20 50 

L=27T/k 

Figure 3. Comparison of power spectra for the three models at a late epoch. This figure shows that non-linear evolution 
has erased large differences between model I (thick line), model II (dashed line) and model III (dot-dashed line). The dotted 
straight line corresponds to a power spectrum with index n = — 1. 
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Figure 4. Growth of power at three scales L = 8, 12, 24 for models I (thick line), model II (dashed line) and model III (dot- 
dashed line). This figure shows that at late times the power in all three models approaches the same value. In particular, the 
growth of power at L = 24 is not influenced by additional power present at smaller scales in model II and III, showing that 
non-linear small scales do not influence larger scales. 

in figure 4. The thick, dashed and dot-dashed hnes represent models I, II and III respectively. 
Curves have been labeled by the length scale. The thick lines demonstrate power generated 
by mode coupling at the two smaller length scales. From the dashed lines for model II, one 
can see that power at L = 8 decreases with respect to the linear rate of growth so as to 
asymptotically match with the amplitude in the reference model, i.e., model I. The same 
effect is seen in the dot-dashed lines for model III. This figure also shows that the existence 
of power at L = 8 does not influence evolution of power at L = 12, if there exists power at 
a larger scale. In this sense, gravitational clustering transfers power from larger to smaller 
scales ["cascades"] but not in the opposite direction ["does not inverse cascade"]. Figure 3 
and 4 demonstrate this fact very clearly. 

3.1 Critical Index 

The two panels of figure 5 illustrate two features related to the existence of fixed points in 
a clear manner. In the top panel we have plotted index of growth Ua = {dln^{a,x)/dlna)x 
as a function of ^ in the quasi-linear regime. Curves correspond to an input spectrum with 
index n = —2, —1,1. The dashed horizontal line at = 2 represents the linear growth rate. 
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An index above the horizontal hne will represent a rate of growth faster than linear growth 
rate and the one below will represent a rate which is slower than the linear rate. It is clear 
that - in the quasi-linear regime - the curve for n = —1 closely follows the linear growth 
while n = —2 grows faster and n = 1 grows slower; so the critical index is Uc ~ —1. The 
curves are based on the fitting formula due to Hamilton et al ( |199lD . Other fitting formulas 
suggested by Jain, Mo and White (|1995| ) and Peacock and Dodds (|1996D give somewhat 
different curves but all these models have fixed points close to = — 1. 

The lower panel of figure 5 shows the slope Ux = —3 — {dlnC,/dlnx)a of C, for different 
power law spectra. It is clear that crowds around rZc ~ —1 in the quasi-linear regime. 

As an aside, we can derive an upper limit on the index of an arbitrary power spectrum. 
Consider the pair conservation equation [eqn.(20) of Nityananda and Padmanabhan (1994)] 
for a set of particles. We can rewrite this equation by using the definition of the index of 
the power spectrum and the index of the rate of growth Ua 

na = h^- (4) 

for the growing mode, the amplitude of perturbations is a monotonically increasing quantity, 
ensuring that the index is always positive. The pair velocity h is also a positive quantity 
at scales with > 0. Therefore the index satisfies the following inequality. 

< I (5) 

We have plotted the curve 3/f in the lower panel of figure 5 as a dotted line. This inequality 
must be satisfied by non-linear structures that have grown out of small inhomogeneities via 
gravitational instability. Any distribution of mass that does not satisfy this inequality could 
not have formed due to gravitational collapse only. 

The index Uc = —1 corresponds to the isothermal profile with ^(a, x) = a^x~^ and has 
some interesting features to recommend it as a candidate for fixed point. For example, in 
the n = —1 spectrum each logarithmic scale contributes the same amount of correlation 
potential energy, [for more details, see Padmanabhan ( [L996a| ).] If the regime is modeled 
by scale invariant radial flows, then the kinetic energy will scale in the same way. It is 
conceivable that flow of power leads to such an equipartition state as a fixed point though it 
is difficult prove such a result in any generality. Equipartition of kinetic energy and the role 
of n = — 1 as the transition index has been pointed out previously by Klypin and Melott 
( |1992|) . They studied evolut ion of kinetic energy at different scales for various models to 
arrive at this conclusion. It would be interesting to see whether the existence of such a fixed 
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a 



a 




Figure 5. The top panel shows exponent of rate of growth of density fluctuations as a function of amplitude. We have plotted 
the rate of growth for three scale invariant spectra n = — 2, — 1, 1. The dashed horizontal line indicates the exponent for linear 
growth. For the range 1 < <5 < 100, the n = — 1 spectrum grows as in linear theory; n < — 1 grows faster and n > ~1 grows 
slower. The lower panel shows the evolution of index rix = —3 — {dln^/dlnx)a with ^. Indices vary from n = —2.5 to n = 4.0 
in steps of 0.5. The tendency for Ux to crowd around ric = —1 is apparent in the quasi-linear regime. The dashed curve is a 
bounding curve for the index {rix < 3/^) if perturbations grow via gravitational instability. 
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Figure 6. Influence of small scales on large scales. The axes are same as figure 2. Thick lines shows the initial power spectrum 
and other curves show the non-linear spectrum at various epochs, labeled by a. This figure shows generation and evolution 
of the fc4 tail. It is seen that the slope of this tail changes rapidly and later goes over to the expected quasi-linear index of 
n = —3/8. Reference lines with these slopes are represented as dashed lines. 

point can be derived, starting from the first principles, say, from the study of the equation 

(S)- 

The eqns. § also show that, in the non-linear regime with ^ > 200, the fixed point is 
'tT'cNL = ~2. Speculating along similar lines, we would expect the gravitational clustering to 
lead to an profile at the non-linear end changing over to x~'^ in the quasi-linear regime. 



3.2 Influence of Small Scales on Large Scales 

Let us now consider the flow of power to larger scales in gravitational clustering. It is well- 
known that the motion of particles conserving momentum leads to a k'^ tail to the power 
spectrum, if the original power was sub-dominant to at small k [( [Zeldovich 1965| ), ( Peebles 



1974 )]. Figure 6 shows the tail for a simulation which had initial power peaked around 
ko = 27t/Lq]Lq = 8. The initial power spectrum was a narrow Gaussian with amplitude 
adjusted so that the peak reaches non-linearity at a = 0.25. In initial stages of evolution, 
there is evidence for a k'^ tail. The amplitude of tail grows as initially, in comparison with 
the linear rate a^, and slows down at later stages. If we assume the quasi-linear evolution is 
governed by equation ^ then the index n = 4 will change to n = —3/8. In figure 6, we have 



© 1996 RAS, MNRAS 000, pHl 



16 J.S.Bagla and T.Padmanahhan 



q 
d 




Figure 7. This figure shows the rate of growth for power at some modes in the fc* tail. For comparison we have plotted a 
dashed line for growth proportional to o? and another for a^/^. This figure clearly shows that at early times the square root 
of the power grows in proportion with a? and at late times the rate of growth drops below a. However, it does not reach the 
expected quasi-linear rate till the end of this simulation. 

plotted a line with this slope for reference. There is some evidence for the slope approaching 
this value; the evolution of slope is definitely in the right direction. [It is clear from lower 
panel of figure 5 that the index of a n = 4 spectrum evolves rapidly, even when ^ < 1. For 
example, the index evolves to less than 3 for ^ = 0.1.] Shandarin and Melott ( |199CI| ) have 
studied the evolution of the tail in detail using two dimensional numerical simulations. 
They also note the change of slope and a slow decline in the rate of growth of the /c^ tail. 

We have plotted the rate of growth for a few modes in the tail as a function of scale 
factor. This is shown in figure 7 that also has two reference lines for growth proportional 
to a? and to a^/^. The latter is the rate of growth in the quasi-linear regime obtained from 
the scaling relation eqn.(^.) It is clear that the rate of growth for all the modes shown here 
starts out as o? but drops to a rate slower than a at late times. 

We stress that the evolution of power outside the band containing initial power is entirely 
due to power transfer by non-linear mode coupling. While at smaller scales this transfer is 
significant and leads to equal amount of kinetic energy per logarithmic wave band, the fiow 
of power to larger scales is less. One can easily see that k~'^/^ spectrum will contribute an 
amount of energy 

^5/8 

per logarithmic band. There is less energy at larger wavelengths, i.e 
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at smaller k. It is, of course, understandable on general grounds that large scales will not 
be affected by the strong non-linearities in the small scales [see e.g., the discussion in §28 of 
Peebles ( |T980D ]. 

4 CONCLUSIONS 

In conclusion, we note that the transfer of power in gravitational clustering shows some 
generic pattern which is worth exploring further. Figure 4 demonstrates that small scales 
[even if highly non-linear] do not influence larger scales. The dominance of cascading over 
inverse cascading as well as the existence of a universal index for the induced the power 
spectrum is reminiscent of fluid turbulence. It may be possible to use some of the concepts 
from the study of turbulence to make ideas like critical indices, fixed points, equipartition 
of energy, etc. sharper and build a new paradigm for understanding non-linear gravitational 
clustering. 
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